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Abstract 

In this paper we present a new semilocal convergence theorem fi'om data at 
one point for the Weierstrass iterative method for the simuhaneous computation 
of polynomial zeros. The main result generalizes and improves all previous ones 
in this area. 
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1 Introduction 

Let / be a monic polynomial of degree n > 2 with simple complex zeros. We consider 
the roots of / as a point in C"". Namely, a point ^ in C" with distinct coordinates is 
said to be a root- vector of / if each of its coordinates is a zero of /. Starting from 
an initial point in C" with distinct coordinates we build in C" the Weierstrass 
iterative sequence [1] 

(1) z^+^ = z^ -W{z^), A; = 0,1,2,..., 

where the operator W in C" is defined by W{z) = {Wi{z), . . . , Wn{z)) with 

W.iz) = /^^j (z = l,2,---,n). 

It is well-known that under some initial conditions the Weierstrass sequence ([I]) is 
well-defined and tends to a root-vector of /. Iteration formula ([1]) defines the famous 
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Weierstrass method (known also as the Durand-Dochev-Kerner-Presic method) for find- 
ing all the zeros of / simultaneously. 

In 1962, DoCHEV [2l[3] proved the first local convergence theorem for the Weierstrass 
method. Since 1980 a number of authors [ai5l|6l[7l[8ll3ll0lllllll2lll3l|lll|15] have 
obtained semilocal convergence theorems for the Weierstrass method from data at one 
point (point estimation). In this note we present a new semilocal convergence theorem 
for the Weierstrass method which improves and generalizes all these results. 

The main result (Theorem [1]) of this note will be proved elsewhere. 



2 Statement of the main result 



Throughout the paper the norm in C" is defined as usual, i.e. \\z\\p = (X^iLi 
For a given point z in C" with distinct coordinates we use the notations 



W{z) 
d{z) 



Wijz) 
di{z) 



Wnjz) 

dJz) 



where d{z) = {di{z), . . . , dn{z)) and 

di{z) = min \ zi — Zj 



and E{z) 



1,2, 



W{z) 



d{z) 



n). 



Theorem 1. Let f he a monic polynomial of degree n > 2 with simple zeros. Let 
1 < p < oo and 1/p + 1/q = 1. Define the real function 



(2) 



{n - ly/'^x 



1 



X 



n.-l 



(1 - x)(l - 2i/9x) V {n- 1)Vp(1 - 21/ax), 
Suppose that z^ is an initial point in C" with distinct coordinates satisfying 
(3) E{z^) < 1/2^/" and (l){E{z^)) < 1. 

Then the following statements hold true. 

(i) The Weierstrass iterative sequence ([1]) is well-defined and convergent to a root- 
vector C, of f . Moreover, the convergence is quadratic if (j){E{z^)) < 1. 

(ii) For each k > 1 we have the following a priori error estimate 



(4) 



e < 



where A = 0(^(z°)) and 9 = 1 - 2^/''E{z°). 
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(iii) For all k > we have the following a posteriori error estimate 

(5) Ik''-' -ell <T^^lk'^'-^1l , 

V ^ II ^iip - 1 _ 0^x1 II lip' 

where Xk = (p{E{z'')) and Ok = 1 - 2^/'iE{z^). 

Remark 1. Let R{n,p) denote the unique solution of the equation = 1 in the 
interval (0,1/2^/'^), where is defined by (|2]). Then the assumption ([3]) of Theorem 1 
can also be written in the form E{z^) < R{n,p). 



3 Comparison with the previous results 



In this section we compare Theorem [T] with all previous results of the same type. Note 
that for z E C"' with distinct coordinates we have the obvious inequality 



(6) 



W{z) 



d{z) 



< 



\\W{z)l 
6{z) 



where S{z) = mm{di{z), . . . ,dn{z)}. 



Corollary 1. Let f be a monic polynomial of degree n > 2 with simple zeros. Let 
1 < p < oo and 1/p + 1/q = 1. Suppose z^ is an initial point in C" satisfying 



(7) 



W{z° 



d{z^ 



< 



1 



A{n - 1)1/9 + 2V9 + 1 



where A = 1.763222. . . is the unique solution of the equation exp {1/x) = x. Then the 
Weierstrass sequence is quadratically convergent to a root-vector of f . Moreover, 
we have the error estimates (jl]) and ([5]). 

Proof. The statement of Corollary [2] holds true even if we replace the right-hand side of 
dZD by = 2/(m + - 46), where m = Aa + b + 1, a = {n - if/i and b = 2^9. 
Define 

9{x) 



ax 



{l-x){l-bx)' 

It is easy to show that g{R) = l/A and 0(x) < g{x) exp g{x) for < x < 1/b. Then 
by the definition of A we get (f){R) < 1 which according to Theorem [1] completes the 
proof. □ 



Batra [13] has proved that the Weierstrass method is convergent under the condi- 
tion II 14^(2°) II oo < S{z^)/ {2n). The following corollary improves Batra's result as well as 
previous results [i El [101 [IH [12] . 
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Corollary 2. Let f he a monic polynomial of degree n > 2 with simple zeros. Let 
1 < p < oo and 1/p + 1/q = 1. Suppose is an initial point in C" satisfying 



W{z°) 



d{z^) 



< 



Then the Weierstrass sequence ([T]) is convergent to a root-vector of f . Moreover, we 
have the error estimates ^ and ([5]). 

Proof. Denote by R the right-hand side of ([8]), i.e. R = l/(2a + 2), where again a = 
(n — 1)1/''. By Theorem [1] it suffices to prove that (j){R) < 1. If n = 2, then (j){R) = 
(4/3) (5 - 6)/(4 - 6)2 < 1, where b = 2^/^ . Further, suppose n > 3. It is easy to see that 



(9) 

If a > 2, then implies 



0(i?)< 



2a{a + 1) 



(2a + l)(2a + 2 -6) 



e. 



2a + 1 5 



If a < 2, then ([9]) implies 



, , 2a(a + 1) r- 3 ^ 

[2a + l)(a + 2) 5 



which completes the proof. 



□ 



Corollary 3. Let f be a monic polynomial of degree n > 2 with simple zeros. Suppose 
that z^ is an initial point in C" satisfying 



W{z^) 



d{z^) 



< R = 0.307541 . 



where R is the unique solution of the equation 



(10) 



X 



1 — X 



exp 



X 



1 — X 



1. 



in the interval (0, 1). Then the conclusion of CorollaryU\ holds forp = 1. 

Proof. It is easy to show that (f){x) < g{x) for < x < 1, where g{x) denotes the 
left-hand side of equation ( fTOl) . Therefore, (f){R) < 1 which according to Theorem [T] 
completes the proof. □ 
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Corollary 4 (Han [T^). Let f be a monic polynomial of degree n >2 with simple zeros. 
Assume 1 < p < oo and 1/p + 1/q = 1. Suppose that 



'111 



d{z^) 



{n - 1)1/9 + 21/9 V {n- 1)1/9 + 21/9 



Then the Weierstrass iterative sequence ([T]) converges to a root-vector of f . 

Proof. Let R denote the right-hand side of fill I) . According to Corollary [2] it suffices 
to prove the inequality R < l/(2a + 2), where a and b are defined as in the proof 
of Corollary [H Write R in the form R = g{T), where g{x) = x{l — ax) and r = 
n(2i/" - l)/(a + b). Since n(2i/" - 1) < 1 we get < g{l/{2a)) = l/(4a) < l/(2a + 2) 
and the proof of the corollary is complete. □ 

Setting p = cxD in Theorem [1] and taking into account ([6]) we obtain the following 
corollary The first part of it is due to Zheng [5j and Petkovic and Herceg [15]. The 
second part of the corollary is due to Zheng |6]. 

Corollary 5 (Zheng [5^, '6] and Petkovic and Herceg [15j). Let f be a monic polynomial 
of degree n> 2 with simple zeros. Let < C < ^ and 



1-C)(1-2C)V 1-2CJ - 6{zO) ■ 

Then the Weierstrass method ([T]) is convergent to a root-vector ^ of f . Moreover, the 
error estimate (111) holds with 6 = 1 — 2C and p = oo. 



Petkovic and Herceg have proved that the Weierstrass method is convergent 
under the condition ||iy(^°)||oo < 5(^°)/(an + 6), where a = 1.76325 and b = 0.8689425. 
Note that for sufficiently large n {n > 13624) this result is an immediate consequence of 
Corollary [H The following corollary improves Petkovic and Herceg's result for all n > 2. 

Corollary 6. Let f be a monic polynomial of degree n > 2 with simple zeros. Under the 
initial condition 



Wiz"^) 



d{z^) 



1 

< 



1.76325n + 0.6869 



the Weierstrass method is convergent to a root-vector ^ of f with the second order of 
convergence. Moreover, the estimates (jl]) and dS]) hold for p = oo. 

Proof. The sequence 0„ = 0(1/(1. 76325n + 0.6869)) is increasing for 2 < n < 132 and 
decreasing for n > 132. Hence 0„ < 0i32 < 1 for n > 2. Now the conclusion follows 
from Theorem [TJ □ 
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Corollary 7 (Wang and Zhao [H]). Let f be a monic polynomial of degree n > 2 with 
simple zeros. Let us assume that 

(12) \\W{z^)\\^<C{n)6{z^) where C(n) = - min(x(l + x)"-^ - 2x). 

Then the Weierstrass iterative sequence converges to a root-vector of f . 

Proof. By Corollary [2] and it suffices to prove that C{n) < l/(2n). This is obvious for 
n = 2 and = 3 since C(2) = 0.25 and C(3) = 0.112. . . Define g{x) = 2x - x(l +a;)"-^ 
It is easy to show that there exists a unique real number t such that < t < 2^/^""^^ — 1 
and g'{t) = 0. One can prove that t < l/{2n — 1) for n > 4. Then by the definitions of 
C(n) and t we get 

C{n) = max g{x) = g{t) = 2{n - 1)^7(1 + 

[0,2l/{"-l)-l] 

which implies C{n) < 1/(3?^) for n > 4. This completes the proof of Corollary [71 □ 

Corollary 8 (Wang and Zhao [8j). Let f be a monic polynomial of degree n > 4 with 
simple zeros. Let us suppose that 

\\W(z^)\\i < C(n) 5(z^) where C(n) = -min ( V""^ - x] . 

x>o \'^^j=i jm J 

Then the Weierstrass iterative sequence converges to a root-vector of f . 

Proof. According to Corollary [3] and ([6]) it suffices to prove that C{n) < 0.3. For a given 
n > 2 define the function /„(x) = Yl^=i TT^f''"''^^ ~ ^' easy to verify that for x > 
we have 

. —, x = {n + l)fn+i{x) - nfn{x). 

This implies /„(x) < fn+i{x). Now by the definition of C{n) we get C{n + 1) < C{n) 
for n>2. Hence C{n) < C{4) = 0.279 . . . for > 4 which completes the proof. □ 

Remark 2. Wang and Zhao [7j have also proved that the Weierstrass method ([I]) is 
convergent under the condition 

(3-2^2). „ 
II 111 n — 1 

Note that for n > 3 this result is an immediate consequence of Corollary [1] and (j6]) . 
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4 Remark on the SOR Weierstrass method 



In this section we are concerned with the successive overrelaxation (SOR) Weierstrass 
method 



(13) 



z>'+^ = z'' -hkW{z''), A; = 0,1, 2,... 



where hk = hk{f) G (0,1] is an acceleration parameter. Petkovic and Kjurkchiev 
[16] have noticed that usually the SOR method (fT3|) is faster if is closer to 1. In 1995, 
Wang and Zhao [8] considered the SOR method ([13]) with defined by 



(14) 



, 0.204378 (5(z'= 
hk = mm < 1, — 



i=l 




According to Corollary [3] we can consider the SOR Weierstrass method (fT3ll with hk 
defined by 



(15) 



hi 



min < 1 



0.307541 



' n 

E 

i=l 



Note that the new hk is closer to 1. Moreover, if our acceleration parameter hk is less 
than 1, then it is greater than Wang-Zheng's parameter hk by more than 50%. 
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